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Peak ground acceleration (PGA) is a critical parameter in 
ground-motion investigations, in particular in  earthquake- 
prone areas such as Iran. In the current study, a new method 
based on particle swarm optimization (PSO) is developed to 
obtain an efficient attenuation relationship for the vertical 
PGA component within the northern Iranian plateau. The 
main purpose of this study is to propose suitable attenuation 
relationships for calculating the PGA for the Alborz, Tabriz 
and Kopet Dag faults in the vertical direction. To this aim, 
the available catalogs of the study area are investigated, and 
finally about 240 earthquake records (with a moment 
magnitude of 4.1 to 6.4) are chosen to develop the model. 
Afterward, the PSO algorithm is used to estimate model 
parameters, i.e., unknown coefficients of the model 
(attenuation relationship). Different statistical criteria 
showed the acceptable performance of the proposed 
relationships in the estimation of vertical PGA components 
in comparison to the previously developed relationships for 
the northern plateau of Iran. Developed attenuation 
relationships in the current study are independent of shear 
wave velocity. This issue is the advantage of proposed 
relationships for utilizing in the situations where there are 
not sufficient shear wave velocity data. 
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1. Introduction 


Several major earthquakes have been occurred in Iranian Plateau that some of them are shown in 
Table 1. Attenuation relationships have an integral role in the analysis of seismic risk [1]. These 
relationships are commonly used to estimate the uncertainty in earthquake motion in 
conventional hazard analysis [2]. Major earthquakes in Kobe (1995), Bam (2003) and Chi Chi 
(1999) have highlighted the importance of the vertical component of the earthquake in the extent 
of damage to diverse building systems [3,4]. Numerous studies have been conducted on the 
importance of the vertical component of earthquake that are mentioned in the following. 


Table 1 
Some major earthquakes in Iran. 
Location (Province) Year Damage (death toll) Local magnitude 
Sarpol-e Zahab (Kermanshah) 2018 over 500 wounded 6.35 
Ahar-Varzaghan (East Azarbaijan) 2012 306 killed & 3037 wounded 6.35 
Mormori(Ilam) 2014 over 250 wounded 6.17 
Taze abad (Kermanshah) 2014 3 killed & 243 wounded 6.17 
Sang Chal (Mazandaran) 1957 1500 killed 6.74 
Buin Zahra (Qazvin) 1962 over 12225 killed & 2800 wounded 6.82 
Dasht-e Bayaz (Sothern Khorasan) 1968 about 10000 killed 6.8 
Korizan Khaf (Khorasan) 1979 420 killed 6.19 
North east of Kish island (Hormozgan) 2022 1 killed & 30 wounded 5:7 
Bakhtar Sarjangal (Sistan and Baluchestan) 2022 0 5.6 


Terzi and Athanatopoulou (2021), by studying a historical building in Greece, concluded that the 
role of the vertical component on the total number of effective building cracks is significant [5]. 
Rahai and Arezoumandi (2008) analyzed a reinforced concrete bridge pier and they showed that 
vertical movements during the earthquake increased the response level and extent of the damage 
to the bridge [6]. They also found that vertical movements result in axial forces acting on the 
bridge columns which in turn, cause instability of hysteresis loops and increase the extent of 
damage. Consequently, in addition to the horizontal component of earthquake, the effects of the 
vertical component of the earthquakes should also be considered in the seismic design of 
structures. Vertical accelerations can change gravity loads that control the in-plane lateral load 
capacity of masonry foundations and affect the out of plane overturning stability of thin-walled 
panels [7]. Noteworthy, higher acceptable performance levels of the steel beams (larger plastic 
rotations) occur when both horizontal and vertical seismic components are considered [8]. The 
majority of researchers believe that it is important to consider the effect of the vertical 
component in high-rise structures and those constructed near faults [9]. 


So far, many efforts have been made by researchers to analyze seismic risk by looking for a 
suitable attenuation relationship between effective parameters of PGA such as magnitude, 
distance, fault depth and geotechnical conditions. Statistical modeling such as regression analysis 
is commonly used to estimate the relationship between two or more variables [10]. Many 
scholars such as Zaré (1999), Ambraseys et al. (2005), Sedaghati and Pezeshk (2017), Soghrat 
and Ziyaeifar (2017) and Ghodrati Amiri et al. (2018) used a variety of regression methods for 
the estimation of PGA [4,11—14]. In addition, optimization techniques such as the genetic 


E. Mohammad Kamareh et al./ Journal of Soft Computing in Civil Engineering 7-3 (2023) 129-142 131 


algorithm (GA) have also been used in the last decades to solve a wide variety of optimization 
problems such as structural damage detection optimization [15-17]. Amit Shiuly et al. (2020) 
used the genetic algorithm and an artificial neural network model to predict the PGA of 
Himalayan region in India [18]. They used the hypocenter distance and magnitude as input 
parameters and PGA records at different distances as output parameters. They compared their 
results to the result of other attenuation relationships and concluded that GA and ANN 
algorithms are capable to predict PGA values. 


Sobhaninejad et al. (2007) used the GA with 586 data records from Europe and the Middle East 
to modify the attenuation relationships proposed by Ambraseys et al. (2005) [12,19]. Also, 
Bagheri et al. (2011) used the GA in Alborz, Central Iran and Zagros then proposed new 
relationships to estimate PGA and peak spectral acceleration (PSA) [15]. Regarding the lack of 
attenuation relationships in the vertical direction in the northern Iran plateau, the aim of this 
study is to extract attenuation relationships in the vertical direction using the particle swarm 
optimization (PSO) algorithm. To this aim, firstly, about 270 earthquake records are used to 
extract the initial model (an explicit equation between PGA and involved parameters), and then 
the PSO algorithm is applied to find the proper coefficients to predict the adequate PGA values 
for northern part of Iran. Finally, the results of developed attenuation relationships in this paper 
are compared to previously developed ones. 


2. Research significance 


In the present research it is tried to develop attenuation models for estimation of vertical PGA 
based on the PSO algorithm for the North of Iran. The importance of using such attenuation 
relationships, is in the seismic design of more resistant buildings and preventing more damage to 
buildings. 


3. Methods 


3.1. Particle swarm optimization (PSO) algorithm 


The PSO algorithm was proposed by Eberhart and Kennedy (1995), which was inspired from the 
food searching behavior of birds (called particles) swarms [20]. The PSO algorithm is a subset of 
the swarm-based algorithms, and swarm-based algorithms are a subclass of nature-inspired 
algorithms, and nature-inspired algorithms are a subset of meta-heuristic algorithms. Meta- 
heuristic algorithms use intelligent strategies to explore a wide search space efficiently. These 
algorithms, in contrast to deterministic algorithms, e.g., the simplex method, because of their 
intelligent searching mechanisms are capable of finding global optimum effectively. In the last 
decades, swarm-based algorithms, especially PSO algorithm, have been widely used in different 
engineering problems [21—24]. The PSO algorithm in comparison to other well-known 
algorithms such as the genetic algorithm (GA) is simpler and easier to apply. Different 
comparative studies showed that the PSO algorithm was superior to other meta-heuristic 
algorithms in finding the final optimum solution/model [25,26]. In the PSO algorithm, 
individuals (birds or particles) in an n-dimensional search space are considered potential 
solutions. The initial swarm/population for each variable in the determined range is randomly 
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generated. In an iterative procedure and by fitness function (error function) the particles are 
evaluated and then the best particle (a particle with the lowest value of error) is introduced 
(called gpest). In addition, each particle keeps its best position (position with the lowest value of 
error) (called ppest) that has been visited by the particle so far. In the PSO algorithm, the velocity 
of each particle that moves towards Ppest and also Gyese changes frequently. If the ke particle in 
an n-dimensional space is considered to be x* = x*,xk, ...,xk, the pyose of this particle will be 
presented by (ppest) x = (Poest)k.1» (Phest)k.2) +» (Ppest)k.n- Lhe best particle among all particles 
is represented by Qpest = Ybest.1»Gbest.2» --» Gbestn [27]. The velocity of the kth particle is 
shown by v* = v*1,v*2,...,v<”". The velocity and position of the particles are set using the 
following equations: 


WEA = wvi + cyrand,. ((Poest)x xf) + c,rand;. (Obese a xf) (1) 
xk, = xk + vk, ,k = 1,2,..., pop (2) 


where pop and t, respectively, indicate the number of particles and iterations, v¥ and v‘,, denote 
the velocity of the k” particle in the ¢” and (t+1)" iterations, respectively. Also, w represents the 
inertia weight coefficient, c; and c2, respectively, show acceleration factors that pull the particle 
towards its Ppest aNd Gpest. Moreover, rand; and rand, respectively, represent uniform random 
numbers in the range of [0, 1] and x¥ and x*,, are the position of k” particle in the t and (t+1) 
iterations, respectively [20,28]. Based on the showed pseudo-code of the PSO algorithm in 
Figure 1, the coding was carried out in MATLAB software R2015a. 


Begin 
1: Objective function, Fy(X), X=(X1, X? ,X3,... X*)" 
2: Initialize each particle's velocity and position randomly, X* (k=1, 2, 3, ..., pop) 
3: Find the p?@5* and g?°5* by objective/error function 
4: while t < MaXjteration do (termination condition) 


5: for k=1:pop 

6 Update particle velocity and particle position by Eqs.1 and 2 
7: Evaluate new solution/location 

8 Check boundary conditions and evaluate new solution/position 
9 Update the p??* and g?°* 


10: end for 

11: end while 

12: Output X* (final solution/model) 
End 


Fig. 1. Pseudo-code of the PSO algorithm. 


3.2. Collected data 


The selected area in this study was the North Iranian Plateau, including Alborz, Azerbaijan and 
Kopet Dag regions, which are parts of the active tectonic regions of Iran (Figure 2a). The 
accelerogram data in the present study were selected from the Building and Housing Research 
Center (BHRC) of Iran, according to the classification of the Standard 2800 (2014), which is 
similar to Eurocode 8 (CEN, 2003) [29,30]. Earthquake records were classified based on the 
shear wave velocity of the soil at a depth of 30 m (V430), i.e., I: Ws30 >750; I: 375<V,30<750; TIT: 
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175<V 305375; IV: V30 <175. Finally, about 240 records with a moment magnitude of 4.1 to 6.4 
and epicentral distances less than 400 km (7 to 359 km) were selected for further analysis. 


Figures 2b and 2c represent the distribution of accelerogram stations and coordinates of the event 
epicenter in the study area. Due to availability of the epicentral distance for all selected data in 
the national accelerogram network, it is used as the distance parameter. About 56.66% of whole 
data have the magnitude of 6-6.4, 29.16% of whole data have the magnitude of 4-5, and 
remained data have the magnitude of 5-6. Noteworthy, in the current study, the proposed catalog 
of Karimiparidari et al. (2014) was used to convert the magnitude data which were not on the 
moment scale [31]. 
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(a) The five active tectonic regions (Mirzaei et al. 1998) [32] 


(b) distribution of accelerogram stations (c) The epicenter coordinates of the selected events 


Fig. 2. Location of the studied region, distribution of accelerogram stations and the epicenter coordinates 
of events. 
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3.3. Error evaluation criteria 


To explore the more efficient models for attenuation relationships, different evaluation criteria 
such as root mean squared error (RMSE), mean absolute percentage error (MAPE), mean error 
(ME), coefficient of determination (R’) and adjusted R? (R? Adjusted), are used to assess the results. 
The RMSE and its combination with MAPE are used as the PSO objective/error function (Fy) for 
developing attenuation relationships. 


Furthermore, the Log-likelihood (LLH) statistical test is used for consistency check between the 
measured and calculated data [33,34]. The LLH method, in order to compare two models, has 
widely been used in engineering sciences in the last few years. The value of the LLH criterion 
changes between 0 and 1, which values toward 1 show the high consistency and vice versa. 


3.4. Formulation of optimization objective function 


The considered objective functions in the current work, are presented by Eqs. 3 and 4. It should 
be noted that RMSE minimization was first considered as the objective function (Eq. 3). 
Regarding the results, the objective function in Eq. 4 was used to improve the optimization 
results where a and f are factors determined by sensitivity analysis. Also, the PGA prediction 
equations are proposed after sensitivity analysis of several different functions based on the least 
squared error rate. 


YL (PGAops —PGAcal)” 
n 


Objective Function: Minimize (Fs) = RMSE = (3) 


Objective Function: Minimize (Fy) = (a x MAPE + B X RMSE) (4) 


where PGA,q; is the peak ground acceleration calculated from the attenuation relationships in 
Eqs. 5 and 6. Eqs. 5 and 6, respectively, estimate PGAs based on the soil and fault types (strike- 
slip and thrust-reverse faults). 


O7 
logy = a, + a2 exp(a3M,,) + ay exp(asR) + dg exp(a7F rp) + dg exp(AgFss) + ayy e)%0 (5) 
S30 


logit) = a, + az exp(agMy) + ay exp(asR) + a6 exp(a7V 30) (6) 


where a, to a, are decision variables obtained from the PSO model, Y stands for PGA, Mw is 
the moment magnitude, R represents the epicentral distance (km), V,, shows the shear wave 
velocity for the top 30 m of soil (m/s), Fs is the strike-slip fault parameter; Fp is the thrust- 
reverse fault parameter and oy,,, and y,,, , respectively, show the standard deviation and mean 
of shear wave velocities of different soil groups at a depth of 30 m. 


4. Discussion 


Sensitivity analysis was performed to determine the PSO parameters. According to the results, 
the internal parameters of the PSO algorithm, i.e., population size, maximum iteration, 
acceleration coefficients (c; and c)) and inertia weight coefficient (w) were set equal to 300, 
1000, 2 and 1, respectively. In this article, different forms of objective functions and attenuation 
relationships were tested. Finally, the best results were selected for presentation. The objective 


E. Mohammad Kamareh et al./ Journal of Soft Computing in Civil Engineering 7-3 (2023) 129-142 135 


function in Eq. 3 was used to obtain attenuation relationships for PGA, estimation. The results 
include the best obtained values and the standard deviation of residuals (Table 2). As seen in 
Table 2, the accuracy of the attenuation relationships is improved by classifying the soil based on 
the shear wave velocity. Table 3 presents the values obtained for the decision variables, i.e., the 
coefficients of attenuation relationships for different soil types with different shear wave 
velocities. 


Table 2 
The optimization results for estimating PGA, in the case of using RMSE as the objective function (Eq. 3) 
Performance measures 


RMSE, MAPEL”  Gvesidua 


The best value obtained for objective function 


eK 


I 0.114 

Soil groups | II 0.327 0.22 79.842% 0.251 
I 0.716 

All input data 0.292 0.292 80.256% 0.355 


"RMSE; is the weighted mean of RMSEs obtained for every group 
“MAPE; is the weighted mean of MAPEs obtained for every group 
oO residual 18 the standard deviation of residuals 


Table 3 
The optimum values of attenuation relation cofficients for estimating PGA, in the case of using RMSE as 
the objective function 


Values of decision variables (coefficients of the attenuation relations) 


ay a a3 a4 as a6 a7 ag ao 410 ait 


I 0.7438 -7.5173 -0.1525 3.1234 -0.0054 -0.2762 -0.0812 -0.0907 -1.9434 1.6421 -0.9170 
Soil groups | IT -0.5187  -9.0524 -0.333 3.0502 -0.0152 0.2737 -1.7204 0.7468 -4.5034 0.8724  -2.9745 
Il 3.2382 -6.6793 -0.2334 1.2775 -0.0118 -3.1572 -3.8436 -2.9659 -5.0928 1.7035 2.4601 
Allinput data -0.8350 -10.00 -0.4256 2.0401 -0.0137 7.0245 -4.3108 0.2116 0.5903  -0.3303 -4.1501 


The general objective function, a xX MAPE + 6B X RMSE (Eq. 4), was used to improve the 
accuracy of attenuation relationships for PGA, estimation. After sensitivity analysis based on the 
different values of coefficients a and £, the objective functions MAPE+2RMSE and 
MAPE+3RMSE were, respectively, selected for the Tabriz-Kopeh Dagh and Alborz-Khazar 
faults. In Table 4 acceptable results are shown. Comparing the results in Table 2 and 4 indicates 
an improvement of MAPE using the combined objective function. Table 5 lists the values that 
obtained for the coefficients of the attenuation relationship for different fault types. 


Table 4 
The optimization results for estimating PGA, in the case of using RMSE and MAPE as the hybrid 
objective function (Eq. 4) 


The best value obtained for objective Performance measures 
function MAPE MAPE,~— ‘RMSE: Oresidual 
Tabriz-Kopeh 0 
Fault types Dagh ue Scag 42.682% 0.316 0.288 
Alborz-Khazar 1.718 46.422% 


"RMSE; is the weighted mean of RMSEs obtained for each group 
™MAPEr is the weighted mean of MAPEs obtained for each group 
6 residual is the standard deviation of residual 
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Table 5 
The optimum values of attenuation relation cofficients for estimating PGA, in the case of using RMSE 
and MAPE as the hybrid objective function. 
Values of decision variables (coefficients of the attenuation relations) 
ay a a a4 as a a7 
Tabriz- Kopeh Dagh 0.19602 -8.7819 -0.2684 2.1508 -0.01614 -4.0111 -9.5617 
Alborz- Khazar __-0.7448 9.9998-  -0.4495 1.4815 -0.01645 2.4233 -3.5297 


Fault Types 


4.1. Comparing the proposed attenuation relationships with other ones 


Derived attenuation relationships in this study were compared with those presented in the 
literature for the north plateau of Iran. To this aim, five attenuation relationships (Table 6) that 
have been proposed for Europe-Middle East and Iran (North Iranian Plateau) were selected to 
compare with the proposed relationship in this study. 


Table 6 
Candidate attenuation relationships. 
Attenuation Models Abbreviation Type of Distance M,, Range 
Zaré (1999) [11] Z-9 Focal 3.2-6.4 
Nowroozi (2005) [35] N-5 Epicentral 5.1-6.4 
Ghodrati Amiri et al. (2007) [36] G-7 Focal 5.1-6.4 
Soghrat and Ziyaeifar (2017) [4] SZ-17 Epicentral 4.1-6.4 
Sedaghati and Pezeshk (2017) [13] SP-17 Boore-Joyner 3.2-6.4 


The evaluation indices including R’, R? Adjusted, ME, RMSE, MAPE and LLH were calculated for 
different attenuation relationships and the results were presented in Table 7. In this table, PGA,-1 
is the proposed model based on the shear wave velocity and objective function of RMSE for 
PGA, and PGA,-3 is the proposed model based on the fault type and hybrid objective function of 
MAPE and RMSE for PGA,. According to the results, PGA,-1 has a lower RMSE than the other 
selected relationships. The PGA,-1 and PGA,-3 models with a higher R’ and R’ adjusted 
outperformed than other relationships. The PGA,-3 and SZ-17 models showed lower MAPE 
values than other ones. However, PGA,-3 outperformed SZ-17 in terms of other evaluation 
criteria (Table 7). Based on ME, the PGA,-3 and PGA,-1 predicted PGA, values lower and 
higher than their recorded values, respectively. Moreover, PGA,-3 with a lower LLH 
outperformed the selected attenuation relationships. The lowest LLH was obtained for PGA,-1. 
The attenuation relationships are ranked in Table 8 based on different assessment criteria. 
According to Table 8 and all assessment criteria, it could be concluded that the proposed 
attenuation relationships in this article outperform other attenuation relationships. 


Table 7 
The values of error evaluation criteria for PGA, estimation relations 
Models MAPE% RMSE ME R° Ro vajisted LLH 
“PGA, -1 (Present Study) 79.843 0.22 0.012 0.593 0.586 0.049 
“PGA, -3(Present Study) 42.682 0.316 0.081 0.999 0.991 0.386 
SZ-17 41.623 0.365 0.073 0.152 0.136 0.612 
SP-17 50.670 0.425 0.117 0.039 0.013 0.815 
N-5 55.441 0.342 0.069 0.209 0.202 0.500 
Z-9 66.267 0.376 0.108 -0.024 -0.043 0.740 
G-7 103.569 0.639 0.078 0.179 0.138 1.400 


‘The proposed model based on the shear wave velocity and objective function of RMSE for PGA, 
“The proposed model based on the fault type and hybrid objective function of MAPE and RMSE for PGA, 
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Table 8 
Ranking of PGA, estimation models. 
Models MAPE% RMSE ME R Radistea LLH 
“PGA, -1 (Present Study) 6 1 1 2 2 1 
“PGA, -3(Present Study) 2 2 5 1 1 2 
SZ-17 i 4 3 5 5 4 
SP-17 3 6 7 6 6 6 
N-5 4 3 2 3 3 3 
Z-9 5 5 6 7 7 5 
G-7 7 7 4 4 4 7 


‘The proposed model based on the shear wave velocity and objective function of RMSE for PGA, 
The proposed model based on the fault type and hybrid objective function of MAPE and RMSE for PGA, 


For further investigation and a better comparison of the proposed attenuation relationships with 
other attenuation relationships, the observed and estimated PGA, values are presented in Figure 
3. As seen, there is greater compliance between the accelerogram and calculated data for the 
PGA,-3 model. 


For further assessment of the attenuation models to each other, they were investigated and 
compared in terms of accuracy and standard deviation of residuals in relation to the moment 
magnitude (Mw). The results are presented in Table 9 where P, and P,, respectively, represent the 
p-values of the slope and y-intercept. The null hypothesis is rejected when the p-values are lower 
than or equal to the significance level of 5%. A p-value greater than 5% leads to confirmation of 
the null hypothesis. When the p-values approach 1, the resulting attenuation model with a lower 
standard deviation will be capable to predict PGA with the higher accuracy. 


cme of residuals versus moment magnitude for PGA, estimation models. 

Models Caner (Pi) Medsedeeiduat 

PGA, -1(Present Study) 0.354 0.410 

PGA, -3(Present Study) 0.743 0.385 

SZ-17 0.575 0.957 

SP-17 0.519 0.244 

N-5 0.162 0.068 

Z-9 0.894 0.422 

G-7 0.258 0.224 


“The P-value of fitted line slope on the residuals (distribution of models residuals versus My) 
The P-value of fitted line intercept on the residuals (distribution of models residuals versus My) 
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(b) The attenuation relation proposed by Zaré (1999) 
[11] 
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(d) The attenuation relation proposed by Sedaghati 
and Pezeshk (2017) [13] 
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(f) The attenuation relation proposed by Ghodrati 
Amiri et al. (2007) [36] 


Fig. 3. The comparison of the attenuation relationships for estimating PGA,,. 
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According to Table 9, Z-9 followed by PGA,-3 has the lowest standard deviation due to a high P- 
value. However, as shown in Table 7, PGA,-3 has a significantly lower ME than Z-9 and thus is 
superior to Z-9. Figure 4 shows the line fitted on the residuals versus the moment magnitude for 
the best and worst PGA, estimation attenuation models. 


PGA,-3 G-7 
5 4.5 
4 
4 7 . 
3.5 
3 a 
= S 
22 
 , — ¥=-0.0085x + 0.1289 ° a 7 
e oe, @ . 
0 Steteeete tit Ht y = -0.1645x + 1.0536 
0.5 ° . 
1 a 
4. AS 5 55 6 6.5 5 5.5 6 6.5 
M, M,, 
a) The best attenuation relation b) The worst attenuation relation 


Fig. 4. Residual distribution versus M,, for the best and the worst PGA, estimation relations. 


Considering MW=6 and Vs30=500 m/s, the PGA, was estimated for different epicentral 
distances using different attenuation relationships (Figure 5). The PGAv-3 output has the 
maximum consistency with the SZ-17 and N-5 output in the vertical direction. There is a small 
gap between the SP-17 diagram and other mentioned diagrams (the type of distance in this 
relation is different from the other relations). The G-7 due to the lack of accelerograph data is 
very different from other charts. 

———S SP 17 = = HN 


mmm a S717 cccccce PGV3 
100 


PGA, (m/s?) 


Epicentral Distance(km) 


Fig. 5. A comparison between the attenuation models for the horizontal and vertical directions considering 
M,=6, V.30=500 m/s (SC=ID) for vertical direction. 


5. Conclusion 


In this study, regarding the importance of the vertical component of PGA in the extent of damage 
to buildings, PGA estimation in the vertical direction was considered. About 240 earthquake data 
were gathered from the BHRC with a moment magnitude of 4.1 to 6.4 and the epicentral distance 
of less than 400 km. Then, using the PSO algorithm, new attenuation relationships were 
developed for seismic regions in the northern plateau area of Iran. Using different evaluation 
criteria and the LLH test, the proposed attenuation relationships were compared with those 
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presented in the literature for the north plateau of Iran, Europe and the Middle East. The obtained 
results are summarized in following: 


PSO algorithm can be used as a powerful optimization method to improve the accuracy of 
attenuation relationships. 

The hybrid objective function considered in this study is very effective in improving the 
accuracy and efficiency of the proposed attenuation models. 

In contrast to the previous attenuation relationships presented by other researchers, some 
attenuation relationships proposed in this study do not need the shear wave velocity data 
for the estimation of ground motion acceleration. This characteristic facilitates their 
application in regions where the shear wave velocity data is not available. 

The attenuation models proposed in this paper outperformed other attenuation relationships 
in terms of different evaluation criteria and they were relatively highly accurate in the 
estimation of the PGA. 

The model proposed for estimating PGA, with the combined objective functions and based 
on the classification of fault types was highly consistent with the model presented by 
Soghrat and Ziyaeifar (2017). 
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